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The possible stable singular vortex (SV) and half-quantum vortex (HQV) of the superfluid 
3 He-A phase confined in restricted geometries are investigated. The associated low-energy exci- 
tations are calculated in connection with the possible existence of Majorana zero modes obeying 
non-Abelian statistics. The energetics between those vortices is carefully examined using the 
standard Ginzburg-Landau (GL) functional with a strong-coupling correction. The Fermi liq- 
uid effect, which is not included in the GL functional, is considered approximately within the 
London approach. This allows us to determine the stability regions in pressure, temperature, 
and applied field for SV and HQV. The existence of the Majorana zero mode and its statistics, 
either Abelian or non-Abelian under braiding of SVs, is studied by solving the Bogoliubov- 
de Gennes equation for spinful chiral p-wave superfluids at sufficiently low temperatures. We 
determined several conditions controllable external parameters for realizing the non-Abelian 
statistics of Majorana zero modes e.g., pressure, field direction, and strength. 

KEYWORDS: chiral p-wave superfluids, vortex, zero energy state, non-Abelian statistics, Ginzburg-Landau 
theory, Bogoliubov-de Gennes equation 



1. Introduction 

Much attention has been paid to various exotic vortices 
and the associated low-energy excitations in spin-triplet 
p-wave superfluids and superconductors. 1-3 > In partic- 
ular, there has been a considerable number of investi- 
gations for the quasiparticle whose energy eigenvalue is 
exactly zero and whose creation and annihilation opera- 
tors are self-Hermitian: 7_e =u = 7g = o- 4 ' 5 '* Fermionic field 
operators that have self-Hermitian creation and annihi- 
lation operators are called the Majorana fermions, 6 -* and 
their quasiparticle energy modes with je^q = Je=o are 
called the Majorana zero-energy states (ZESs). The Ma- 
jorana ZESs localized at vortices are considered useful 
for fault-tolerant quantum computations because they 
can obey the non-Abelian statistics. 5, 7 ~ 10 ) 

The localized Majorana ZESs have been pointed out in 
the quasi-hole of the Pfafnan state of the quantum Hall 
state with the 5/2 filling. 4,5 -* The other candidate sys- 
tems that support Majorana ZESs are quantized vortices 
in chiral p-wave superfluids or superconductors, 7, n ~ 17 ) 
the surface Andreev bound state of chiral or time re- 
versal symmetric p-wave superfluids, 18 ~ 23 ) the junction 
between an s-wave superconductor and a topological in- 
sulator, 24 ~ 27 ) and the s-wave superfluid with particular 
spin-orbit interactions. 28 ^ 

First, we consider Majorana ZESs bound at a singu- 
lar vortex in spin-polarized p-wavc superfluids, which 
involves the particle- hole symmetry je = 1-e- Then 
the existence of the exact ZES is guaranteed by the 
index theorem 29 ) and the analytic solution of the 
Caroli-de Gennes-Matricon state when the vorticity is 
odd. 12,17,30,31 -' Therefore, the existence of the Majorana 
ZESs is topologically protected. The eigenstates of the 
system constructed by Majorana quasiparticles are de- 
scribed by the occupation of the complex fermion state, 
that is, a linear combination of two Majorana quasipar- 
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tides, c 2 j = (72.7-1 +ij2j)/y/2. For spin-polarized sys- 
tems, a complex fermion is formed by a pair of the Ma- 
jorana ZESs "fj and Jj+i in between spatially separated 
vortices, giving rise to the nonlocality of the complex 
fermion. If a well-isolated vortex has an odd number of 
Majorana ZESs, the complex fermion state is necessar- 
ily constructed by the spatially separated Majorana ZES 
and the vortices obey the non-Abelian statistics because 
of the nonlocality of their excitations. 7 - ) For the system 
where the spin degrees of freedom, for example, the spin 
Iff) and I44-) pairs, survive, an isolated singular vortex 
(SV) has a degeneracy of the ZESs in two spin sectors. 32 ) 
The statistics of the SV in spinful superfluids depends 
on the ways of forming a complex fermion as follows. 
If the two Majorana quasiparticles constructing a com- 
plex fermion eigenstate arc in a single vortex core, the 
statistics is Abelian. If they are in spatially separated 
vortices, however, the statistics is non-Abelian. Thus, a 
half-quantum vortex (HQV), that has a singularity in 
either the spin Iff) or \H) component of the order pa- 
rameter, is more feasible for the non-Abelian statistics. 7 ^ 
The superfluid 3 Hc-A phase is one of the most effective 
candidates of the p-wavc superfluid associated with the 
Majorana ZES in the sense that the ground states are 
most established 1,33 ' 1 and there are many experiments 
and samples with controllable parameters. 34 ~ 38 ) Since 
the superfluid 3 He-A phase is a spinful system, the statis- 
tics of the vortex in spinful p-wave superfluids should be 
clarified. 

In order to observe the Majorana character, it is neces- 
sary to set up a two-dimensional system where the thick- 
ness of the sample is much less than the dipolc coherence 
length. A possible experiment has been set up by Saun- 
ders group, where they confine superfluid 3 He in a slab 
geometry whose thickness is 0.6 /zm under the follow- 
ing conditions: The pressures are between and 5.5 bar, 
the temperature is cooled to 350 /iK, and the external 
field is applied perpendicular to the slab to carry out the 
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NMR observation. Then the group observes the phase 
transition from the A phase to an unknown phase in a 
low-pressure and low-temperature region. Investigation 
of vortices in spinful chiral p-wave supcrfluids and their 
features will also lead us to understand such a character- 
istic of 3 He-A in a slab in future. 

For rotating experiments, Yamashita et aZ. 36 ) have re- 
cently performed an experiment on parallel-plate geom- 
etry intended to observe the HQV in superfluid 3 Hc-A. 
The superfluid is confined in a cylindrical region with a 
radius R = 1.5 mm and a height of 12.5 /xm, sandwiched 
between parallel plates. A magnetic held H = 26.7 mT 
(||z) is applied perpendicular to the parallel plates, and 
the pressure P=3.05 MPa. At this pressure, the strong- 
coupling effect due to spin fluctuations becomes impor- 
tant. 39 > Since the gap of 12.5 /J,m between plates is small 
compared with the dipole coherence length £<j ~ 10/xm, 
the Z-vectors, which indicate the direction of the orbital 
angular momentum of Cooper pairs, are always perpen- 
dicular to the plates. Also the d-vectors are confined 
within the plane by the applied held H \\ z because the 
dipole magnetic held is Hd ~ 2.0 mT, 1 ) where H tends to 
align the d-vectors perpendicular to the held direction. 
Yamashita et al. conducted an investigation in various 
parameter spaces, such as the temperature T, or the ro- 
tation speed fl up to f2 = 6.28 rad/s using the rotating 
cryostat in ISSP, University of Tokyo, capable of achiev- 
ing a maximum rotation speed of ~T2 rad/s; however, 
there is as yet no evidence of the HQV. 36 ) 

The stability of HQVs has been argued since HQVs 
were pointed out by Volovik and Mineev. 40 ) The hydro- 
dynamical calculation taking account of the Fermi liq- 
uid (FL) correction shows that the HQV is energetically 
stable against the SV 2,41 ~ 44 ) However, in our previous 
study, 45 ) we have carried out the calculation on the basis 
of the full Ginzburg-Landau (GL) theory taking account 
of the strong-coupling correction due to spin fluctuations, 
which play a crucial role in the stability of the A phase 
at high pressures, without the FL correction. This im- 
plies that the HQV is unstable unless the pairing phase 
becomes the so-called A2 phase under a strong external 
held. 45 ) In this study, we hnd that the contributions of 
the strong-coupling and FL effect to the stability of the 
HQV are competitive at realistic rotation speed and high 
pressure in the experiment using the rotating cryostat in 
ISSP. 

Therefore, in order to clarify the statistics of vortices 
in the spinful p-wave superfluid, we have to examine the 
statistics of the SV, in addition to that of the HQV on 
an equal footing. The argument on the statistics of vor- 
tices is valid in adiabatic and quantum limits, where the 
ZESs are energetically distinct from the other core bound 
states with finite energies. This implies the absence of the 
decoherence of ZESs and requires the temperature to be 
sufficiently lower than the energy difference between core 
bound states T^/Tp where T c is the transition tempera- 
ture and Tp is the Fermi temperature. 

This paper is arranged as follows: In §2, we derive the 
spinful Bogoliubov-de Gennes (BdG) equation and intro- 
duce the GL free energy. In §3, we discuss the energetics 
of the vortex textures and the stability of the HQV. We 



calculate the energetic advantage of the HQV originating 
from the FL correction and the disadvantage originating 
from the spin-fluctuation strong-coupling correction us- 
ing both the London limit calculation and the full GL 
calculation. In §4, we examine the statistics of the SV 
using numerical calculation with the spinful BdG equa- 
tion. In §4.1, we discuss the case in which the direction 
of d- vectors is perpendicular to the external field. In this 
section, we clarify that when the vortex distance is fi- 
nite, two eigenstates originating from ZESs in different 
spin sectors do not hybridize with each other so that the 
braiding of the SV does not commute necessarily. In §4.2, 
we consider the case where the direction of d-vectors is 
tilted from the direction perpendicular to the external 
field. In this situation, we demonstrate that the Zeeman 
effect due to an external field parallel to the d-vectors 
causes the hybridization of the two Majorana ZESs. In 
the hnal section, we present our summary and conclu- 
sions. 

2. Formulation 

2.1 Spinful Bogoliubov-de Gennes equation 

In general the OPs for spin triplet superfluids are de- 
scribed as 



A(n,r 2 ) 



A tt (n,r 2 ) A n (r 1 ,r 2 ) 
A4.-1- (n,r 2 ) A u (n,r 2 ) 



(1) 



Htl'i,'-2j 

where A aa > = -V(r\, r 2 ) (V^CrOVv (r 2 )) and ip a (r) is 
the held operator of fermions with spin a =f, 4- The spin- 
ful mean-field Hamiltonian of the spin-triplet superfluids 
and superconductors is described using this notation as 



H = E + ^Jdr 1 J dr 2 * t (r 1 )£(ri J ra)*(r 2 ) 



¥(r) = f^t(r),^(r),^}(r),^{(i 



(2) 



(3) 



H (ri,r 2 ) A(ri,r 2 ) 
-A*(n,ra) -flj(ri,ra) 



, (4) 



H {ri,r2)=S(r 1 -r 2 )[H {r 1 )l + ti n H-d-}, (5) 

where ^(r) is the spinor in the Nambu space, 1 is a 2 x 2 
unit matrix, and /Lt n , H, and <r are the magnetic moment 
of 3 He atoms, the external held, and the 2x2 Pauli ma- 
trices respectively. Here, we assume that the superfluids 
and superconductors are confined by the potential with 
a magnitude Vq 3> lA^'l. Then the single particle part 
Ho(r) is given as 



Ho(r) = - 



2M 



V 9(r -R)-fi + ifl-(rxV), (6) 



where M, /x, and fi are the mass of the particle, the chem- 
ical potential, and the external rotation, respectively. We 
set H = 1. 

In general, the OP of p- wave superfluids is expanded in 
terms of the eigenstate of the orbital angular momentum 
of the Cooper pair l z = —1,0,1. We assume the system to 
be spinful chiral p-wavc supcrfluids, namely, the orbital 
ferromagnetic state, so that the only l z — — 1 component 
survives. This can be realized in a parallel-plate geometry 
and a slab, where the dipole coherence length is much 
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larger than the thickness of the sample. Then the explicit 
expression of the OPs is given as 



n + n 



, k I — A c 



ri + r 2 



ki 



■exp[-(fc 2 -fc 2 )£ 2 ], 



(7) 



where k is the relative momentum and £ p is the size of 
the Cooper pair. The components in eq. (1) are obtained 
as the Fourier transformation of eq. (7) with respect to 
the relative coordinate 



ixi2 + Vl2 
x — : — — exp 



4£ 



(8) 



where kp is the fermi wave number. 

We carry out the Bogoliubov transformation from the 
fcrmionic field operator to the quasiparticle basis with 
r„ = [IVf-, r„4,, T^, T^] T in the Nambu space: 



(2) r (i)\* r (2)V 



where r i is the 4x4 matrix defined by a Pauli matrix. The 
transformation matrix u v (r) must satisfy the orthonor- 
mality J dr\ul(ri)u^(ri) = cV, p l and the completeness 

J2u *"iM„( r l)^4( r 2) = S(ri - r 2 )l, where 1 = diag(i, 1) 
is a 4 x 4 unit matrix. In order to obtain the operators of 
the quasiparticles that diagonalize the mean field Hamil- 
tonian eq. (2) as 

h = e + frt^iv, (9) 



(10) 



one can find that the wave functions of the quasiparticles 
ti^(r) and u^\r) should obey the same BdG equation 
described as 

J dr%(r, r')u v (r') = E^u v (r), (11) 

u v {r) = [ul{r), ui(r), vl(r), «£(r)] . (12) 

Then the annihilation operator of the quasiparticle is ex- 
pressed as 



dr {w l/ (r)} t *(r). 



(13) 



We numerically diagonalize the BdG equation (11) un- 
der the gap potential given in eq. (7) with 



4^,-i(r) = ^ctct'.-i II ex PKA'] tanh ( J > ( 14 ) 



i=i 



where Kj is the winding number of the component in 
the j-th vortex, 9j and rj are the azimuthal angle and 
radius centered by jth vortex, respectively, and £ aa i 
is the coherence length of the OP in the spin sec- 
tor \<jcr'). We assume the uniformity of the OPs along 
the z-direction. For the quasiparticle eigenstate in the 



BdG equation (11), we impose the periodic boundary 
condition with the wave number k z . that is, u v = 
UE,k z (x, y) exp[ik z z] and T va = TE,k z ,a- Then the BdG 
equation (11) is block-diagonalized in terms of k z . In the 
subspacc, the particle-hole symmetry {z x IC z x }* = — £ 
gives Z x { u E,k z exp[ik z z]}* = u_p.k z exp[i/s z z] and the 
inversion symmetry along the z-direction gives up^ r = 
up,-k z and E(—k z ) = E(k z ). Thus, one finds Tp k a = 
T-s.-fc-.a, which implies that the quasiparticle arising 
from k z = can be the Majorana zero mode rj a = 
Po,o,<t = 7 CT - We focus on the eigenstate with k z = 
throughout this work. The numerical diagonalization 
is carried out by the discrete variable representation 
method. 46 - 49 ' 

2.2 Ginzburg- Landau framework 

We use the GL framework to discuss the stable tex- 
tures and the energetics of the vortex, which is quanti- 
tatively reliable for 3 He. Then we assume that the OPs 
are decomposed to the center-of-mass coordinate and the 
orbital degrees of freedom with a relative momentum 
around |fc| ~ kp. These components are described with 
the 3x3 matrix = A^i{r) as 



A(r,fc) 



-A x i + iAyi \plA Z i 
\2<A Z i A x i + iAy 



A^m 
A-\lrn 



A^^rn 
Allm 



fei(15) 



(16) 



where i = x,y,z, m = — 1,0, +1, k± = (k x =p ik y )y/2, 
and k is the unit vector oriented to the direction of the 
momentum on the Fermi surface. The GL free-energy 
functional, which is invariant under gauge transforma- 
tion and under spin and orbital space rotation, is well- 
established 1 ' 2, 33 ~ 35 > 5 °- 53 ) an d given by the standard form 

/total = /bulk /bulk /grad + /dipolc + /field; (17) 



/bulk — ~ a ^* f ii^iy 



/bulk = PiA^A^AvjAvj + l3 2 A* tll A* uj A^ i A 1J j 
Jr fi-j,A* lli A* vi A jl jA V j + ^A^AljA^j Aui 
iA* ,A vi A v 



(18) 



(19) 



/ g rad = KtftA^WArt) + K2{d*A* M ){d 3 A^) 

+K 3 (d*A; i )(d j A„ j ), (20) 
/dipoic = Sd(^^i- + A^Ay^ — T^A^A^), (21) 

/field = 9mH ^A^HvAvi. (22) 

In the weak-coupling limit, the GL parameters and cou- 
pling constant of the dipole energy 52 ' are 

a = a (l - T/T c ), a = AM 



0? = f>Y = 0? 



-PY = -20? 



_ 0R W _ 7C(3)jV(0) 
— Z 'I J ~ 120(7rfc B T e ) 2 ' 

K _ K _ K _ 7C(3)jV(0)(to f ) 2 
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3d 



40 



7ftAr(0)lnLi 339xp2> 



7C(3)JV(0)( 7 ft) 2 
™ 48[(l + F °)irfc B T c J • 

The details of the physical constants in the form de- 
scribed above are as follows: the transition temperature 
T c , the density of states N(0), the Fermi velocity vp, 
the permeability of vacuum /xo, the gyromagnetic ra- 
tio 7, and the Fermi temperature Tp. These are given 
by the experiments 34 ' 35 ' and depend on the pressure. In 
the high-pressure region, the GL parameters of the bulk 
4th-order terms ft are corrected by the strong-coupling 
effect due to spin fluctuations. For ft, we use the strong- 
coupling correction calculated by Sauls and Serene, 50 - 1 as 
mentioned below. 

3. Energetics of Vortex Textures 

Without the loss of generality, we use the description 
for the OP = d^Ai, where the d- vector d^ and Aj 
are complex values. For the bulk of the 3 Hc-A phase, d- 
vectors can be unit vectors, but we consider the generic 
situation on d-vectors throughout this work. First, let us 
define textures under the situation that all the d- vectors 
lie in the xy-plane. This situation is approximately real- 
ized in 3 He between parallel plates under the strong field 
H ^> Hd- If we choose the direction perpendicular to the 
plane as the spin quantization axis, then A^ jTn = 0. In 
this case, we can find the following two possibilities on 
the vortex textures. One of them is the HQV: In the vor- 
tex core, either the A^ jm or Ay. >m component of the OP 
has a unit winding number and singularity. The other is 
the SV: The singularities of both the spin components of 
OPs are in the same position. 

In the zero-field case, a single SV has two zero energy 
modes originating from the spin and |++) sectors, 
whereas the HQV has a single zero energy mode from 
the spin sector that has a phase singularity in the OP. 
The low-energy bound state at the vortex core of the 
HQV is the same as that of the singular vortex in spinless 
chiral p-wave superfluids. Therefore, it is well-known that 
HQVs obey the non-Abelian statistics. 7 ' 

However, the energetics of the HQV against the SV 
still remains a problem. Specifically, for the realistic set 
of the GL parameters in 3 He, the FL correction maintains 
the stability of the HQV. In contrast, it is demonstrated 
here that the strong-coupling effect on the bulk 4th-order 
terms that stabilizes the A phase against the B-phase is 
not favorable for the stability of the HQV. On the basis 
of the GL theory, we examine here the stability of the 
HQV as a consequence of the competition between the 
strong-coupling effect and the FL correction. 

The strong-coupling effect in the bulk 4th-order term 
is derived by Anderson and Brinkman 39 ' as 

ft = —(I + 0.1<5)ftT, ft = (2 + 0.2<5)ftf , 

ft = (2 - 0.05<5)ftf, ft = (2 - 0.555)ftf , (23) 

ft = -(2 + 0.75)ftf, 

where 6 is the spin-fluctuation parameter depending on 
the pressure. 33 ' Thus, the bulk 4th-order term of the GL 



free energy (19) is 

/bl = £ d (Mtt| 4 + \ d u\ A ) + B c (\d n \ 2 \d U \ 2 ), (24) 
where 

B d = ftf [(4 - 0.355)(|A +1 | 4 + \A^) 

+ (16-0.55<5)|A+ 1 | 2 |A_ 1 | 2 ] , (25) 



B c = -P™5 [3.5(|,4 + i | 



l^4-i| 4 ) 



+9L4+ 1 | 2 |A_ 1 | 2 ] <0, 



(26) 



where A m = (A x - sgn(m)iA y )/V2, d aa = (d x - 
sgn(cr)idj / )/v2, and g^j, = d z . If both the spin |ft) an d 
\H) components of the OP remain finite, the free energy 
decreases with increasing 5 because B c < in eq. (24). 
The amplitude of the OP in both the spin sectors is en- 
hanced by each other through the strong-coupling effect. 
In the case of the core of the HQV, one of the spin com- 
ponents must have the singularity so that another com- 
ponent is not enhanced by the mechanism arising from 
B c < 0. 

In order to quantify the strong-coupling effect, we 
carry out the numerical minimization of the GL free en- 
ergy in eq. (17) composed of the bulk, gradient, and 
dipole energy in eqs. (18)-(21). Our numerical condi- 
tion is as follows. We assume that the I- vectors align 
to the z-direction and the d-vectors are in xy-plane. 
The uniformity of the OPs is also assumed along the 
z-direction so that the spin and orbital indices of OPs 
Api reduce to fx, i = x, y and our calculation is carried 
out in a 2D plane. Then the magnetic interaction en- 
ergy in eq. (22) can be ignored. This situation can be 
realized in a parallel-plate geometry, 36, 37 ' where the dis- 
tance between the parallel plates is less than the dipole 
coherence length (~ 10 /an) and much greater than the 
coherence length (~ 10 nm). Here, we impose the rigid 
boundary condition with the radius R on the OPs as 
^^(M = R) = 0- We use the GL parameters ft in 
eq. (19), taking account of the strong-coupling correc- 
tion given by Sauls and Serene, 50 ' and their values are 
qualitatively consistent with eq. (23). In Fig. 1, we show 
the amplitudes of the dominant components A-^^i and 
for the SV and HQV. In order to compare these 
energies on an equal footing, we set up two singular vor- 
tices for the SV texture and four half-quantum vortices 
for the HQV texture. In our calculation, however, the 
HQV texture is only the saddle-point solution of the GL 
equation since the HQV can be continuously transformed 
into the SV texture and the free energy of the HQV is 
always higher than that of the SV as long as the FL 
correction is neglected. 

In refs. 42-44, the OPs in the HQV are restricted 
within the bulk A phase. The strong-coupling effect 
equally affects the free energy of both the HQV and SV 
texture under this assumption. In our work, we inves- 
tigate the energetics of the vortices without restrictions 
even at the vortex core and take account of the strong- 
coupling effect. Thus, the strong-coupling effect near the 
vortex core gives the HQV relative energetic disadvan- 
tage compared with the SV, as shown below. 
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Fig. 1. (Color online) Spatial profiles of the dominant component 
|-<4<t<t,+i| for the four-HQV (a) and the two-SV (b) textures near 
the center of the system, and (c) their cross section magnifying 
the vortex core region. The upper (lower) panel in (c) is the spa- 
tial profile of the HQV (SV). The nonsingular component of the 
HQV is depressed at the vortex core signified by the arrow in (a) 
and (c). The broken line in the upper panel of (c) is the ampli- 
tude of |A CTIT +i| of the SV when the position of the vortex core 
coincides with that of the HQV. In all the figures, the system 
size is set to be R = 5 /im, the rotating speed is Q = 3 X 10 3 
rad/s, and temperature T = 0.95T C , where the SV is energeti- 
cally stable. The unit of the x- and y-axis is micrometers and 
the amplitudes of the OPs are normalized with 7rfcsT c . 



Wc analyze the energetics between these two textures 
in detail. In Fig. 2, we show the difference of the bulk 
and gradient energies of the HQV from those of the SV, 
normalized with |A| function of the pressure. 

The features of the HQV texture result from two fac- 
tors: (i) As seen in Fig. 1(c), the amplitude of the com- 
ponent j4^|_ j+ i is depressed at the singularity of A^+i. 
(ii) Near the singularity, the component A^+i in the 
HQV is more enhanced than that of the SV, as shown 
in the upper panel of Fig. 1(c), because the bulk 4th- 
order term with B c < in eq. (26) makes the interaction 
between A^ >m and An im attractive. The bulk energy 
of the |LJ.) component is increased by the factor (i) and 
that of the spin |ft) component is reduced by owing to 
the factor (ii). As seen in Fig. 2, the bulk energies of the 
HQV relative to the SV is negative in the high-pressure 
region, and as the pressure decreases, that is, the strong- 
coupling correction becomes less important, the bulk en- 




20 

Pressure [Barl 



Fig. 2. (Color online) Difference of the gradient and bulk energies 
of the HQV from that of the SV normalized with ii"|A| 2 as a 
function of the pressure. 



ergy of the HQV eventually becomes higher than that 
of SV. This result implies that the bulk energy with the 
strong-coupling correction favors the HQV. 

However, the gradient energy of the HQV is consider- 
ably larger than that of the SV at vortex cores. As seen 
in Fig. 2, the gradient energy of the HQV relative to that 
of the SV becomes larger than that of the bulk energy, 
so that the HQV is not stable with this consideration 
without the FL correction. The reason for this disadvan- 
tage of the gradient energy for the stability of the HQV 
is as follows. By the factor (i), that is, the depression 
of the OP without the vortex singularity at the vortex 
core of the HQV seen in Fig. 1, the loss of the energy 
appears owing to the spatial modulation. Furthermore, 
the kinetic energy due to the phase winding is also en- 
hanced by the fact that the singular component A-^-^i 
is enhanced by the factor (ii) . The strong-coupling effect 
plays a crucial role in the stability of the HQV in the 
high-pressure regime, while this becomes less important 
in the low-pressure regime so that these disadvantages 
due to the gradient energy may be negligible. 

Although the FL correction is the key factor for the 
stability of the HQV, 42-44 ' it is difficult to take account 
of these effects into our calculation on the basis of the GL 
framework using a more generic form of OPs. In order to 
discuss the FL correction, we apply the London approx- 
imation to the OPs, where A + i(r) = \Aa\ cxp[i<f>(r)], 
A-i(r)=0, and the e?-vcctor is assumed as the unit vec- 
tor [d x ,d y ,d z ] = [cosa(r),sina(r),0]. We can treat the 
FL correction as the effective mass of the spin current 
using this representation. Note that this representation 
restricts the pairing phase to the A phase so that we can- 
not describe the vortex core and boundary of the sys- 
tem exactly. The FL correction makes the mass of the 
spin current less than that of the mass current, that is, 
Psp/Ps < 1 where p sp and p s are the effective masses of 
the spin and mass currents, respectively. 

The gradient term with the FL correction is described 
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by the London approximation as 33 ' ' 

m^^r m+rxn)2+ipMiVa)2 - (27) 

If we consider an axially symmetric vortex, $ = <j s </>, 
a = q sp 4>, where <p is the azimuthal angle of the system 
centered in the vortex singularity and q s and q sp are the 
mass and spin circulations, respectively. The HQV tex- 
ture has the circulation (q s ,q sp ) = (1/2,1/2), whereas 
the SV has the circulation (<f,<f p ) = (1,0). Thus, the 
spatial variation term of the GL free energy for each tex- 
ture is 

f dr 

f K 1 + ^) ln (&) ( HQV ) 

= \ (28) 

where the energy of the HQV texture is twice larger than 
that of the single axially symmetric HQV texture, and 
we set the external rotation to be = 0. For p sp = p s , 
the energy of the HQV arising from the spatial varia- 
tion of the OPs is equivalent to that of the SV. If p sp is 
smaller than p s owing to the FL correction, the energetic 
advantage of the HQV texture increases logarithmically 
as the system size R increases. 

We numerically estimate this effect in the London limit 
that neglects the vortex core under the same conditions 
as those for our full GL calculation, such as the geome- 
try and vortex configuration. From the GL calculation, 
we find that the amplitude of OPs recovers to the bulk 
within about 0.5 (im from the vortex core, as shown in 
Fig. 1. Then we neglect this region from the contribution 
and calculate the spatial variation energy using cq. (27). 
Our numerical calculation is carried out as follows. The 
phase factor is taken as the axially symmetric form for 
the vortex core $ = q s J2i 4>i> a = Q sp J2i where fa is 
the azimuthal angle for the j-th vortex core. The geom- 
etry of the system and the layout of the vortices are the 
same as the GL calculation shown in Fig. 1. We set the 
rotating speed to be il = ft c where the two-SV texture 
becomes energetically stable against the one-SV texture 
at the system size R. For example, when the system size 
R = 100 fim, the critical speed is found to be il c = 7.5 
rad/s which is the feasible rotating speed in experiments 
using a rotating cryostat in ISSP. 36 ' The integral on the 
left-hand side of eq. (28) is carried out numerically, and 
the distance of the vortices is determined by minimizing 
the gradient energy. This distance is consistent with the 
calculation based on the full GL theory. Our model is 
so simple that it is sufficient for our purpose to estimate 
only the energy scale of the advantage of the HQV. 

As seen in Fig. 3, the HQV advantage is on the or- 
der of 10~ 32 , normalized with K\A\ 2 p s / p®. By compar- 
ing Figs. 2 and 3, we find that the relative instability 
of the HQV originating from the strong-coupling effect 
becomes comparable to the energetic advantage of the 
HQV when the FL correction is taken into account. 

The strong-coupling correction in the bulk 4th-order 




Fig. 3. (Color online) The main panel is the difference of the 
spatial variation energy / gra d normalized with K\A\ 2 p B /p® as 
a function of the system size R where the FL correction is set to 
be psp/ps = 0.9. The rotating speed is set to be f2 = Q c where 
the two-SV texture becomes stable compared with the one-SV 
texture for p Bp /p B = 1. In our calculation, f2 c is determined as 
a function of the system size shown in the bottom inset. The 
energetic stability of the HQV is proportional to 1 — p Bp /p B as 
seen in the top inset. 



terms decreases as the pressure decreases, and the FL ef- 
fect becomes important as the temperature decreases. 
Hence, the HQV is energetically stable in the lower- 
pressure and lower-temperature region. However, in or- 
der to carry out the quantitative calculation of the en- 
ergetics, we need the GL formulation systematically in- 
cluding the strong-coupling bulk 4th-ordcr terms and the 
FL correction of the gradient terms on an equal footing 
with the general form of the OPs. If we use the most 
general representation of the OPs, that is, the spin and 
orbital parts of OPs are not separated, it is impossible 
to introduce the FL correction by the phenomenological 
method used by Cross. 41 ) These quantitative calculations 
of the stability remain as a future problem. 

As shown above, the HQV has a single ZES bound at 
the vortex core and obeys the non-Abelian statistics. On 
the other hand, SVs are energetically comparable to the 
HQV. Thus, we consider their statistics in the following 
section. 

4. Excitations and Braiding of Singular Vortex 

As shown in the previous section, the stability of the 
HQV remains a problem. Thus, in this section, we con- 
sider the structure of the excitation and braiding of the 
SV. In particular, we notice the Zeeman effect due to 
an external field, and consider the following two situa- 
tions: The magnetic fields applied perpendicular to the 
d-vectors and tilted from its direction. In this section, 
we clarify how the spin degrees of freedom of ZESs af- 
fect the statistics of vortices. Here, we only focus on the 
ZESs and ignore the contributions of the other low-lying 
core-bound states, because the coupling with the other 
core-bound states due to thermal fluctuations may give 
rise to the decoherence of ZESs. This requires the system 
that we consider here to be in the quantum limit, where 
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the temperature is lower than the energy difference be- 
tween core-bound states, that is, T^/Tp. For 3 He, since 
T c 2 /Tp ~ 10 -6 K, the temperature regime assumed in this 
section corresponds to the sub-micro Kelvin. This regime 
does not correspond to the temperature region that we 
considered within the GL framework in the previous sec- 
tion. However, we consider this limit as a starting point 
to discuss the statistics of vortices in realistic 3 He sys- 
tems. In addition, in the case of ultracold atomic gases 
with a p-wave Feshbach resonance and polar molecules, 
the superfluid transition temperature T c is close to Tp, 
because the pair interaction is flexible under an external 
field. In this case, the braiding operation and the real- 
ization of the statistics of vortices are more feasible. 

4-1 H perpendicular to d-vectors 

The excitation composed of the self-Hcrmitian opera- 
tors j2j-i an d 72j describes the complex fermion eigen- 
statc defined as 



C2j = (723-1 + H2j)/V2. 



(29) 



The many-body ground state is described by the oc- 
cupation number of the complex fermion. The complex 
fermion state can be understood as a spatially nonlocal 
state when the self-Hermitian excitations 723—1 and 723 
are localized at spatially separated vortices. If the en- 
ergy eigenvalue of the excitation is exactly zero, an adi- 
abatic braiding of a vortex around another one changes 
the occupation number. 7 ~ 9 ) For instance, let us consider 
that there are four vortices Vi (i = 1,2,3,4) and they 
have one Majorana quasiparticle 74 at each vortex. We 
assume that the Majorana quasiparticlcs 71 and 72 (73 
and 74) form the complex fermion c 2 (4). If the vortex 
V 2 moves around V3 , the operator of this braiding is de- 
scribed as 7,8 ^ 



T23 



(c| + c 4 )(4 - c 2 ), 



(30) 



where we ignore the phase factor. This operator changes 
the occupation number of the complex fermion c 2 j. This 
feature is known as the non-Abelian statistics of the vor- 
tices. When the odd number of the ZES appears at each 
vortex core, not less than one complex fermion are com- 
posed of the ZESs between different vortices. Then these 
vortices obey the non-Abelian statistics. 

On the other hand, we consider the case where the 
number of ZESs localized at one vortex core is even. 
For instance, the ZESs are spin-degenerate at the core 
of the well-isolated SV discussed in the previous sec- 
tion. This leads to the even-number degeneracy of the 
self-Hermitian operators 7J and 7^ at the core of each 
SV Vi, arising from the spin Iff) and \H) sectors. We 
consider that there are four vortices. One can find two 
possible ways of forming a complex fermion. (i) The self- 
Hermitian operators at the same vortex can make the 
complex fermion described as 



Oi = (rfl + i^)/y/2. 



(31) 



In this case, an exchange of the vortex means an exchange 
of the complex Dirac fermion Oj, implying that the braid- 
ing operator gets the only phase factor -1 and belongs to 
the Abelian group, (ii) The self-Hermitian operators at 



the spatially different vortices make the complex fermion 
described as 



b2i = il2i 



H 2 , 



(32) 



In this case, if the vortex V2 moves around V3, the braid- 
ing operator is described as 



T23 - ({b\} f + bl)({blY-bl) 

x({bi} f +bi)({biy-bi) 



(33) 



When this transformation operates the quasiparticle vac- 
uum state, the four complex fermions b\, b\, b\, and b\ 
are created. Although the realization of b^ is not pro- 
tected topologically, the braiding operator can change 
the occupation of the complex fermion and the vortices 
obey the non-Abelian statistics. 

The situation where the ZESs are exactly degenerate 
can have both ways of forming the complex fermions at 
and &2r However, the degeneracy of the ZES is removed 
by the finite distance of the vortices through quasipar- 
ticle tunneling. We find that the complex fermion state 
is constructed from the ZESs belonging to different vor- 
tices in one of the spin sectors when the external field 
is oriented exactly perpendicular to the direction of the 
d-vectors as shown in this section. We will show that a 
self-Hermitian particle belonging to a spin state does not 
hybridize with its counterpart in different spin sectors. 
Therefore, the braiding of the vortices changes the occu- 
pation number of a complex fermion in each spin sector 
and the braiding operator is approximately non-Abelian. 

Under the assumption that the external field H is ap- 
plied to the direction perpendicular to the d-vector, we 
can block-diagonalize the BdG equation (11) into the two 
spin sectors |ff) and \\\) as 



dr 2 



where, 



K. n (r l7 r 2 ) 

t n {r u r 2 ) 



Uu v {r 2 ) 
= EJAu v {r 1 ), (34) 



-A* CT (r 1 ,r 2 ) -F -(n,r 2 ) 



Uu v {r) = [ul,vl,ut,v^ 



(35) 
(36) 



where U_ is an appropriate 4x4 unitary matrix and the 
spin quantization axis is set to be the direction parallel 
to H . Notice that eq. (35) is equivalent to the Hamilto- 
nian density of spinlcss p-wave supcrfluids and the com- 
ponents of the wave function belonging to |ff (||)) arc 
given as (vt™, t/lw). In the case of the SV, the two 
ZESs appear in two spin sectors, namely, |ff) and |44)- 
We consider systems with a plural number of SVs. For 
an infinite vortex distance D v , the ZES originating from 
each sector degenerates precisely. Any linear combination 
of these two states can be the eigenstate of the system. 
However, for a finite D v , this degeneracy of ZESs can be 
removed by the following two factors. One is the Zeeman 
effect for energy splitting due to quasiparticle tunneling 
and the other is the splitting of the coherence length of 
two spin components under a strong external field. 
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First, we discuss the former factor. In a spinless p- 
wave supcrfluid, two ZESs bound in neighboring vortices 
tunnel and interfere with each other. Then the energy 
of the complex fermions constructed in an intervortex 
increases from zero so that the complex fcrmion states 
are not exactly ZESs. 49 ' 54 ) This energy shift oscillates 
and decreases exponentially as D v increases, originating 
from the quantum oscillation and localization of the wave 
function of ZESs. 

We apply previous spinless argument in the spinful 
case. When we take account of the Zeeman effect, the 
chemical potentials of the up-spin and down-spin parti- 
cles effectively shift as /i±/j, n H, where \i n is the magnetic 
moment of the particle. The energy shift of ZESs in the 
spin |ff ) a nd sectors can be calculated separately as 

2|^° 



E„ 



r 3/2 



cos(fc L> v + V 4 ) 



± \/D v ko sin(Z? v k 

k 



Cp 



o 



tt/4) 



exp 



£,cro 



(37) 



where fen = 2Mfx — Because of the second term in 
this expression, E a , the eigenvalues of the spin Iff) sector 
and |44) sector deviate from each other, and they can- 
not hybridize with each other. Thus, complex fermions 
should be constructed between vortices in each spin sec- 
tor. 

In order to demonstrate this, we numerically diago- 
nalize the spinful BdG equation (11). This diagonaliza- 
tion is carried out in a two-dimensional system under 
the OP given in eqs. (8) and (14). The system has a cir- 
cle geometry and the supcrfluid is confined by the rigid 
wall potential, as in the GL calculation. Here, we set 
|^ct CT _iI > O.Icf in eq. (14) in order to secure the dis- 
creteness of the eigenvalue within the accuracy of calcu- 
lation. Here, we consider two vortices that are located 
at (x, y) = (-Dy/2,0) and (D v /2,0). Although this cal- 
culation is carried out in the strong-coupling region in 
the sense of \A® a _ 1 | ~ ep, the features of ZESs are in- 
dependent of the details of the Hamiltonian. Therefore, 
the result of this calculation can be qualitatively applied 
to the weak-coupling supcrfluid such as 3 He. 

In Fig. 4, we show the wave function of the first- and 
second- lowest excitations in the presence of Zeeman split- 
ting, fi n H/eF = 1.0 x 10~ 4 , at a distance of the vor- 
tices kpD v = 20. These excitations are the complex 
fermion state as a consequence of the tunneling of the 
ZES bound at vortex cores. Then near a vortex core, 
the self-Hermitian relation u° = {i>^}* appears approx- 
imately. As shown in Fig. 4, these excitations are ap- 
proximately composed of only one spin component of u° 



and v°. The u\ and v\ components of the first-lowest 
excitation are on the order of 10 -7 of the components 
u\ and v\. This means that the eigenstates originating 
from the |ff) and |44-) sectors are well-separated. We 
present in Fig. 5 the lowest eigenenergies, which oscillate 
and decay exponentially as a function of D v . In Fig. 5, 
the difference between the eigenvalues of the two diffcr- 
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(Color online) Cross sections of the amplitude of the wave 
Mj, and vi at y = 0. We set the system size 
kpR = 30, the Zeeman splitting fj, n H/tp = 10~ 4 , the coherence 
length of the spin \cr<j) component kp^a- = 2.5, and the distance 
of the vortices kpD v = 20. The panel in (a) shows the wave 
function of the lowest excitation with Ei/ep = 2.08 X 10 — 4 and 
the panel in (b) shows that of the second-lowest excitation with 
E 2 /e F = 2.09 X 10" 4 . 




Fig. 5. (Color online) Lowest eigenenergies E\ and E2 originating 
from the spin |"|~|~) and |4_|.) sectors and the difference E\ — E2 
as a function of the distance of vortices. Here, we set the system 
size kpR = 30, the Zeeman splitting [i n H/ep = 10~ 4 , and the 
coherence length of the spin \crcr) component kp£ arT = 2.5. The 
left and right axes indicate the scale of the eigenenergies E a and 
the difference E\ — E2 , respectively. 



ent spin sectors is found to be about 10 -2 times larger 
than the amplitude of eigenvalues E a . and to have the 
phase difference 7r/2 from the eigenvalue oscillation, as 
shown in cq. (37). In realistic cases, since the amplitude 
of the OP is much smaller than the Fermi energy and 
the distance of the vortex is larger than this numerical 
simulation, these separations of the eigenvalue are quite 
small but finite. Therefore, even in a realistic situation, 
we conclude that complex fermions can be constructed 
approximately by Majorana ZESs between different vor- 
tices in each spin sector. 

In the case of spinful chiral p-wave supcrfluids under a 
strong external field, the amplitude and coherence length 
of the OP component of the spin |ff ) pair are not equal 
to those of the \\X) pair. In the language of 3 Hc, this situ- 
ation is called the A2 phase. The coherence length t; aa of 
the dominant spin component becomes smaller than that 
of the minor component and the ZES is tightly bound 
at the vortex core. The interference through tunneling, 
shown in eq. (37), is weak for the major component and 
strong for the minor one. Therefore, the energy difference 
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Fig. 6. (Color online) Lowest eigenenergies E\ and E2 originating 
from the spin (solid line) and |44-) (dashed line) sectors as 
a function of the splitting rate X defined in eq. (38), where we 
set k F R = 30, fj, n H/e F = 1CT 4 , and k F D v = 16, 18, 19, and 
20. The energy differences of the eigenstatc originating from two 
spin sectors (solid and broken line) are on the order of 10 — 6 €f 
at the splitting rate X = 0. They are enhanced by increasing |X| 



to the order of 10 ep- 



between the ZESs of each spin sector is enhanced by the 
deviation of from The eigenvalue in the domi- 
nant spin sector is closer to zero than that in the minor 
one. Thus, the eigenstate in the former has the Majo- 
rana character more precisely in the sense that the ex- 
citation of the complex fermion state is degenerate with 
the vacuum state. The excitation spectrum continuously 
approaches to that of the spinlcss case as this imbalance 
of OPs increases. 

We numerically diagonalizc the BdG equation for the 
OP given as 

|Aa| 



l^n.-il 



1^2 



u, 



V2 
V2 



il + X) 



(1-X) 



(38) 



In a realistic system, the splitting X 6 [— 1, 1] is pro- 
portional to the external field. 55 -* In Fig. 6, we show the 
behavior of the energy splitting as a function of the im- 
balance of the OP. The solid lines are the energy of the 
eigenstate in the sector E 2 and the dashed lines are 
that of the spin |ft) sector. In the region X > (X < 0), 
the eigenvalue E\ (E 2 ) in the spin |ft) sector be- 

comes larger and the other eigenvalues E 2 (Ei) become 
closer to zero so that the difference E\ — E 2 also increases. 
For instance, at a splitting rate X = 0.48, this difference 
is on the order of 10 _3 £f- In addition, this energy split- 
ting is more enhanced by the short distance of vortices, as 
shown in Fig. 6. The component closer to the zero energy 
has a Majorana character more precisely because of the 
degeneracy with the vacuum state of the quasiparticle. 

It is not desired that the energy difference between 
the occupied and vacuum states of the complex fermion 
is finite. Because of the energy difference, the adia- 
batic exchange of vortices cannot change the eigenstate 
of the complex fermions. Thus, we should carry out 
the braiding for the finite time scale eF/\A° aa _ 1 | 2 < 



t < (ii^-) -1 , where the lower boundary depends on the 
energy discreteness of the core bound states, AE ~ 
|^4ctct -1 | 2 / e F- 17 ' 1 In the realistic value of 3 He, the time 
scale of (So-) -1 is almost infinite and ep/lA^ _ 1 \ 2 is also 
much larger than the time scale of experiments. In the 
case of a p-wave resonant Fermi gag 11 - 13 . 15 - 17 ) r two- 
dimensional polar fermionic molecules, 56 ~ 58 ) this range of 
the time scale may be feasible in experiments. Further- 
more, since there is an energy difference in the many- 
body ground state between the states before and after 
the braiding operation, the thermal relaxation changes 
the eigenstate transformed by the braiding. 

4-2 H tilted from H _L d-vector 

In the previous section, we assume that the (i-vectors 
are perpendicular to the external field H . However, it is 
difficult to precisely align the d-vectors to this direction 
in experiments. Thus, we consider the situation where 
the (i-vectors are tilted from the direction perpendicular 
to H. 

Here, we also assume that the (i-vectors are spa- 
tially uniform. Thus, we choose the directions of the 
(i-vectors and the magnetization as d = (1,0,0) and 
H = (H. x , 0, H z ), respectively, where one finds Aff ttn = 

~ A-ll^m- 

We first consider the eigenstate arising from a well- 
separated SV and regard the Zccman effect arising from 
H x as perturbation. In the unperturbed case, the BdG 
equation is block-diagonalized as shown in eq. (34), and 
the wave functions of the ZESs in these two sectors are 
described as 

1 T 





tti(r) = 
u 2 (r) 



uftr), 0, {ul(r)Y 
0, 4(r), 0, {4(r)Y 



(39) 
(40) 



<(r) = exp(i&aWM*%r) exp(-r/£™)- (41) 

Here, & a is the phase of OP A aa ,-\ at the vortex core 
that arises from the phase of the other vortices, Jq(x) is 
the Besscl function, J\f is the normalization constant, and 

= J 'k F + + sgn(a)/j, n H z /eF. Note that fcg" ~ ftp 
for kF^aa S> 1 and /i„i? 2 <C £f, and the phase factor 
$ t - $1 = 7r since A n . m = -A^ >m . According to the 
ordinary perturbation theory, the perturbation H x re- 
moves degeneracy as 



1 ( 



m 2 ) + 0(fi n H x /e F ), 



(42) 



V2 



(til - m 2 ) + 0(fi n H x /e F )- 



This wave function means that the operator of the eigen- 
state under a finite H x is described by a t in cq. (31). 
Therefore, the two ZESs originating from different spin 
sectors hybridize in the same vortex core so that the hy- 
bridized state behaves as the Dirac fermion. As shown in 
the previous section, the vortices that have such a struc- 
ture of the excitation cannot induce the non-Abelian 
transformation of the many-body ground state by the 
braiding of the vortices. 
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Fig. 7. (Color online) (a) Schematic diagram of the parallel-plate 
geometry and the definitions of 6jj, 9^, and H x . (b) Angles be- 
tween the d- vector and the external field H as a function of Oh 
calculated in eq. (44) under the absolute values of the external 
fields H = 27, 50, and 100 mT. 



The first-order energy shift due to perturbation is es- 
timated as 



A£« 



dr [u' 1 (r)y (inH^u'^r) 



(43) 

where the 4x4 Pauli matrix = diagfdi, — <jj\. Here, we 
use the asymptotic form of the Bessel function J (x) ~ 
yj2j (irx) cos(x — 7r/4) and assume that (^ct^f)" 1 <C 1 

and Co-<r(fco _ k o) — £,<j<jk F H n H z / e F < 1 using the phys- 
ical parameters for 3 He. In the case of the SV in 3 He 
in the parallel-plate geometry defined in Fig. 7(a), the 
angle of the d- vectors is determined as 



1 

- — tan 
2 



l26> 



H 



s20 



H 



(44) 



.{H d /Hf- 

where Od and Oh arc defined in Fig. 7(a). Equa- 
tion (44) is derived by minimizing the dipole energy in 
eq. (21) and the magnetic interaction energy in eq. (22) 
within the London approximation, where A + i(r) = 
\Aa\ exp[i$(r)], A_i(r)=0, and the d-vectors are as- 
sumed as the unit vectors. The resulting angle Od — On 
becomes maximum at Oh = 7r/4, as shown in Fig. 7(b). 
For instance, using Oh = 7r/4 and H = 27 mT, we esti- 
mate H x = 7 AO x 10~ 2 mT. Then, we find AE^ = 
1.68 x 10~ 8 ep, which implies that the shift is much 
smaller than the gap \A aa \ ~ lO~ 3 e F in 3 He. However, 
the symmetry separating two spin sectors is broken by 
this perturbation, and the complex fermion at in eq. (31) 
is constructed in the same vortex core even for an in- 
finitesimal field. 

Although the perturbation H x hybridizes the excita- 
tions in spin sectors, the interference through the tunnel- 
ing of the quasiparticlc separates the spin sectors shown 
in the previous section. In order to clarify this conflict, we 
diagonalize the BdG equation (11) where Od — On ^ tt/2. 
In Fig. 8(a), we plot the maximum amplitudes of the 
wave functions |m|(t)| and |uf(r)| of the lowest-energy 
eigenstate as a function of Od — On- As shown in the previ- 
ous section, the component |ui(r)| is finite and u\(r) = 
at Od — On = 7r/2, implying that this eigenstate orig- 
inates from the spin |ft) sector. As seen in Fig. 8(a), 
when Od — Oh deviates from 7r/2, the minor component 
|wj(r)| grows rapidly, implying that the two sectors hy- 
bridize with each other. In fact, in the case of the co- 



herence length kF^aa = 1.5, the two components of the 
wave functions become equal at Od — On = (89/180)7r: 
|ui(t")| ~ |u|(r)|. This is consistent with eq. (42). We 
carry out the calculation under various fJ, n H and kp^ aa 
values. Then, we find that the magnitude of the Zee- 
man shift fi n H does not change the behavior of the hy- 
bridization. As shown in Fig. 8(a), the hybridization for 
tilting Od — Oh weakens with increasing coherence length 
kF^aa- These results imply that the energy splitting at 
Od — Oh = normalized with the amplitude of Zeeman 
splitting \E\ — £'2 1/ {UnH) determines the growth rate of 
the hybridization with respect to Od — Oh- In order to 
quantify this, we define the initial slope of the hybridiza- 
tion for the f (4_)-dominant mode described as 



ta) 



d(0 d - H ) 



q 



til) 



d{0 d -0 H ) \\vl {i) \ 




(45) 



where |u^j ax | and I^MaxI are the maximum values of 
the lowest-energy wave functions |uj(r)| and |uf(r)|, 
respectively. In Fig. 8(b), we plot the initial slope at 
Od — Oh = 7r /2 as a function of \Ei — E 2 \/ (fi n H). As 
shown in Fig. 8(b), all the results under different exter- 
nal fields, and the layouts of vortices (one-SV and two-SV 
case) are on the same function, so that we ensure that 
the initial slope of the hybridization depends only on 
this ratio and yields a power law behavior in the region 
\E 1 -E 2 \/(fx n H) < 10- 1 . 

In addition, as shown in the previous section, the split- 
ting of the coherence length, ^ 7^ £44, involves the en- 
hancement of the energy difference of the ZESs. The ini- 
tial slope decreases owing to this splitting and the re- 
sult is presented in Fig. 8(b) in the small splitting region 
\Ei - E 2 \/{pL n H) < 1 corresponding to X < 1CT 2 . With 
increasing \X\ corresponding to \Ei — E 2 \l{[i n li) > 1, 
the initial slope deviates from the line in Fig. 8(b) 

because the initial slope in eq. (45) is not a proper indi- 
cator of the behavior of the hybridization. This is due to 
the spatial expansion of the wave function in the minor 
spin sector. However in the region \Ei — E^j (fi n H) > 1, 
we find S„ < 10 _1 . Thus, we can control the d- vectors to 
be perpendicular to the magnetic field enough to ignore 
the hybridization. 

One can find that the accuracy of the direction of the 
magnetic field and the weak coupling of minor compo- 
nent of the pair potentials A^^ m or A^ m are required 
by the hybridization of two spin sectors. For example, 
for 3 He using the rotating cryostat under a high external 
field, 36 ) the feasible rotating speed is Q ~ 10 rad/s, the 
vortex distance is D v = 50 /im, and the magnetic field is 
H = 10 T. In this situation, with the tilting angle of the 
external field Oh ~ 0.10 deg, we estimate Od — On = 10~ 9 
deg. It is concluded that even if the coherence length of 
the minor component of the OP is Co-o-^f > 10 4 , we con- 
trol the hybridization within |iij(r)|/|iij(r)| ~ 10~ 2 . 
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Fig. 8. (Color online) (a) Maximum amplitude |itjj ax | of the wave 
functions uf (r) as a function of the angle between the d- vector 
and the magnetic field. The dashed line is |«i{ ax | and the solid 
line is |u^[ ax |- (b) The absolute values of the initial slopes STf and 
of the hybridization of the minor spin components \u\\ (open 

symbol) and (filled symbol), which are defined in eq. (45), 
as a function of \Ei — E%\/ (iM n H). 



5. Conclusions 

We have studied the vortices and low-energy excita- 
tions of the spinful chiral p-wave superfluid on the basis 
of the phenomenological Ginzburg-Landau (GL) theory 
and the microscopic Bogoliubov-de Gennes theory. We 
focus on the 3 He-A phase between parallel plates under 
a magnetic field. 

In spinful chiral superfluids, possible candidates of the 
vortex texture are the singular vortex (SV) and the half- 
quantum vortex (HQV). In §3, we have discussed the 
energetics of these textures. The free energy of the HQV 
is decreased by the Fermi liquid (FL) correction. 42 ~ 44 ) 
However, we find that the strong-coupling correction due 
to spin fluctuations affects the energetics of the vortex 
texture at the vortex core. Hence, this correction makes 
the HQV unstable compared with the SV, as a result 
of our full GL calculation, which is not included in the 
discussion in refs. 42-44. We calculate the contributions 
of the FL effect using the London approximation and of 
the strong-coupling effect using the full GL framework 
separately. It has been demonstrated that the latter effect 
becomes comparable with the former under a rotation 



speed of 5-10 rad/s and near the transition temperature 
T c , which is possible for the experiment for 3 Hc using a 
rotating cryostat in ISSP. The quantitative calculation 
taking account of two effects described above has not 
been established and remains a future problem. 

In §4, we have investigated the low-energy excitations 
and the statistics of the SV in the quantum limit, where 
the discreteness of the levels in the vortex core is suffi- 
ciently larger than the temperature. For 3 He, since the 
energy scale of the discreteness is T^/Tp ~ 10 _6 K, it is 
difficult to realize the quantum limit in practical exper- 
iments. Nevertheless, in this section, we have discussed 
the statistics of vortices in spinful p-wave supcrfluids and 
achieved the conclusions summarized below. This may 
be the starting point for further study of the statistics of 
vortices in realistic situations that show the decoherence 
of ZESs and the dissipation of vortex motion. In addi- 
tion, our results are applicable to vortices in p-wave res- 
onant Fermi gases 11-13 ' 15 ~ 17 ) and polar molecules. 56 ~ 58 ) 
It is well-known that the HQV has the Majorana quasi- 
particle ZESs, which are topologically protected. In SVs, 
the spin |ft) and |||) components of the order parame- 
ter (OP) have a singularity at the same position so that 
the SV texture has two zero-energy excitations localized 
at vortex cores. These behaviors change as the angle be- 
tween the d- vectors and the external field H varies. 

In §4.1, we have considered the case where the d- 
vectors are exactly perpendicular to H. In the case of 
3 He between parallel plates, this situation is realized by 
applying a sufficiently strong external field and suffi- 
ciently controlling its direction to be perpendicular to 
the plates accurately. In this situation, the SV has the 
two degenerate ZESs that originate from the spin |tt) 
and I44-) sectors when the vortex distance is infinite. This 
leads to the Abclian statistics of vortices. However, for 
the finite vortex distance, the ZESs split through the in- 
terference of their wave functions. Then, the degeneracy 
is removed and the spin and \H) sectors cannot hy- 
bridize with each other. Thus, the excitation structures 
are found to be the same as those in the spinless case, and 
the SV approximately obeys the non-Abelian statistics. 
In addition, when the amplitudes of the OP components 
split under a high external field, that is, the coherence 
length of the spin |ft) P a h is not identical to that of 
the spin |44-) pair, the splitting of the eigenenergy is fur- 
ther enhanced by the interference of the zero-energy wave 
functions. 

In §4.2, we have considered the case where the d- 
vectors are tilted from the direction perpendicular to H . 
We have found that when the d-vectors are tilted, the 
eigenstates originating from the spin \\\) and \{\) sec- 
tors hybridize intensely and form complex fermions in 
one of the vortex cores. Then our numerical calculation 
demonstrates that the intensity of the hybridization is 
determined by the energy splitting of the eigenstates of 
two spin sectors when d-vector _L H . In order to control 
this hybridization in experiments, one has to carry out 
the experiment near the spin-polarized state called the 
A 2 -phase. 

Finally, in the case of 3 He, we have discussed that the 
statistics of the vortices in spinful chiral p-wave super- 
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fluids with the vortex distance D v depends on external 
parameters as follows. In the low-temperature and low- 
pressure region, the HQV is stable, and the statistics of 
the vortices is non-Abelian. In the high-temperature and 
high-pressure region, the SV is stable as a consequence 
of the strong-coupling effect. When the coherence length 
of either the spin |ft) ° r component of OPs is much 
smaller than that on the order of 10 _1 D V , we cannot 
control the hybridization of the spin sectors. Hence, it is 
found that the statistics of the SVs is Abelian. However, 
for coherence lengths larger than 10~ 1 D V , the statistics 
of the SVs is found to be non-Abelian. 
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